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(54) Method of automatically selecting Oligonucleotide Hybridization probes 



(57) A method is provided for automatically select- 
ing oligonucleotide hybridization probes for detecting a 
mutation causing a mismatch in a DNA duplex. This 
method comprising the steps of: defining the mutation 
to be detected by a DNA sequence of a mutant DNA 
strand and a DNA sequence of a wild-type DNA strand 
corresponding to the mutant DNA strand; defining an in- 
teger variable L1 for a first length of the oligonucleotide 
hybridization probes from the place of mutation in a first 
direction; defining a range from ml to n1 of possible val- 
ues of the variable L1; defining an integer variable 12 
for a second length of the oligonucleotide hybridization 
probes from the place of mutation in a second direction 
opposite to the first direction: defining a range from m2 
to n2 of possible values of the variable L2; selecting one 
or more kind(s) of the oligonucleotide hybridization 
probes from probes which have an one-to-one hybridi- 
zation fit to the mutant DNA strand, probes which have 
an one-to-one hybridization fit to the wild-type DNA 
strand corresponding to the mutant DNA strand, probes 
which have an one-to-one hybridization fit to an opposite 
sense mutant DNA strand which is the complementary 
DNA strand to the mutant DNA strand, and probes which 



have an one-to-one hybridization fit to an opposite 
sense wild-type DNA strand which is the complementa- 
ry DNA strand to the wild-type DNA strand correspond- 
ing to the mutant DNA strand; determining all possible 
hybridization probes from the possible values of the var- 
iables L1 and L2, and the kind(s) of probes; defining a 
thermodynamic nearest-neighbour model for calculat- 
ing a melting point of a certain hybridization probe hy- 
bridized to a DNA strand, the melting point being a tem- 
perature at which a predetermined percentage of a mul- 
titude of identical pairs of said certain hybridization 
probe and said DNA strand is in an annealed state; cal- 
culating, on basis of the thermodynamic nearest-neigh- 
bour model, for all possible hybridization probes a first 
melting point of the respective probe hybridized with its 
complementary mutant DNA strand or opposite sense 
mutant DNA strand, respectively, a second melting point 
of the respective probe hybridized with its complemen- 
tary wild-type DNA strand or opposite sense wild-type 
DNA strand, respectively, and a temperature difference 
between the second melting point and the first melting 
point; and ranking the possible hybridization probes with 
regard to the temperature difference. 
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[0009] This object is achieved by a method according to claim 1 

ES!! ^ Vantage ° US embodiments of method are outlined in the subclaims 2 to 9 
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Figure 1 is a flux diagram of a first embodiment of the new method; 
Figure 2 is a flux diagram of a second embodiment of the new method; 
Figure 3 is a flux diagram of a third embodiment of the new method; 

Figure 4 shows an example DNA dupiex of a 18mer oligonucleotide; there are 17 n-n doublets and initiation en- 
s ergies for termination with G-C and A-T must be considered; the oligonucleotide is non-self-complementary; 

Figure 5 is a diagram of experimentally observed melting points versus melting points calculated with the n-n 
model; the regression was calculated according to Passing and Bablok [33], while the standard error of the estimate 
(Sytx) was calculated with the principal component method; and 

Figure 6 is a contour plot showing the interdependence of the melting temperature shift of 1 9mer probes destabi- 
10 feed by a single mutation versus their GC content; the different frequencies are shown in shadings. 

[0014] As indicated in Figure 1 , a method of automatically selecting oligonucleotide hybridization probes for detecting 
a mutation causing a mismatch In a DNA duplex starts with the step 1 of defining the mutation to be detected by a DNA 
sequence of a mutant DNA strand and a DNA sequence of a wild-type DNA strand corresponding to the mutant DNA 

15 strand. Here, it is unimportant whether the mutant DNA strand and the wild-type DNA strand are sense or antisense 
strands as the sense strand always defines the antisense strand and vice versa. In the following as well as in the claims 
a term opposite sense strand is used to refer to the strand having the sense opposite to the strand dined in step 1 . 
Therefore, opposite sense strand is no other formulation for antisense strand but may also refer to a sense strand. 
[00151 In step 2 of the method integer variables L1 and 12 are defined for a first and a second lengths of the oligo- 

20 nucleotide hybridization probes from the place of mutation in a first and a second direction opposite to the first direction. 
Further, ranges from ml ton1 and from m2 to n2 of possible values of the variables L1 and 12, respectively, are defined. 
[0016] Alternatively, step 2 can include defining a first variable for an end point of the oligonucleotide hybridization 
probes being located from the place of mutation in a first direction; and defining a second variable for a length of the 
hybriddization probes in a second direction opposite to the first direction. Then, first and second ranges are to be 

25 defined for the first and the second variables, in which the second range is dependent of the first range as the hybrid- 
ization probes have in any case to cover the place of mutation. 

[0017] At least in the end, both alternatives of step 2 have the same results. Therefore, they are in the following as 
well as in the claims regarded as identical. 

[0018] In step 3 one or more kind(s) of the oligonucleotide hybridization probes are selected from 

30 

probes which have an one-to-one hybridization fit to the mutant DNA strand, 

probes which have an one-to-one hybridization fit to the wild-type DNA strand corresponding to the mutant DNA 
strand, 

probes which have an one-to-one hybridization fit to an opposite sense mutant DNA strand which is the comple- 
35 mentary DNA strand to the mutant DNA strand, and 

probes which have an one-to-one hybridization fit to an opposite sense wild-type DNA strand which is the com- 
plementary DNA strand to the wild-type DNA strand corresponding to the mutant DNA strand. 

AH of the above kinds can be selected, if no special requirements are given. If such requirements exist, special kinds 
40 of the oligonucleotide hybridization probes can be selectively selected or deselected. 

[0019] In step 4 all possible hybridization probes are determined from the input data of steps 2 and 3, i.e. from the 
possible values of the variables L1 and L2, and the kind(s) of probes. 

[0020] In step 5 a thermodynamic nearest-neighbor (n-n) model for calculating a melting point of a certain hybridi- 
zation probe hybridized to a DNA strand is defined. That n-n model will be discussed in detail below. The melting point 
45 is a temperature at which a predetermined percentage of a multitude of identical pairs of said certain hybridization 
probe and said DNA strand is in an annealed state, i.e. hybridized. 

[0021] In step 6, on basis of the thermodynamic nearest-neighbor model, the following values are calculated for ail 
possible hybridization probes: 

so a first melting point of the respective probe hybridized with its complementary mutant DNA strand or opposite 

sense mutant DNA strand, respectively, 

a second melting point of the respective probe hybridized with its complementary wild4ype DNA strand or opposite 
sense wild-type DNA strand, respectively, and 

a temperature difference between the second melting point and the first melting point. 

55 

[0022] At last, in step 7, the possible hybridization probes are ranked with regard to the temperature difference, and 
a ranked list of the hybridization probes is output 

[0023] According to Figure 2, a preferred embodiment of the new method further includes step 8 in which a temper- 
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[0024] As the maximum Tm shift caused by the mutation is limited ftisalsosufffcwt,.,^ 

at least on Tphysico^hemL Seter^ ^S>Z£ ^ Paramet6r ° f ,he "Ration assay. The 
which the hybniaton oHhVhyZSon ?J£X5Z? TIT * ™ ,0n ^"ration under 
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neither able to hybridize to *— »**- ""i** are 
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^ri 1^1 8 bac ^— d ^en for the n-n Mode, to be defined «n step 5 of tt^e methods according to 

Kratu^en^ 



T m=—J- 273.15 



AS [/VfH-]+/?xln(CT) ^ 1 1 

" m is the enlropy at a g ,en 
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oligonucleotide: 

= I^H°(/) + AH°(51nit) + AH°(3W) + AH°(symmetry) [Eq 2] 



AH 0 , 



4 



EP1 103 910 A1 



AH 0 (4 is the enthalpy change for the ten possible n-n and n, is the number of occurrences of the respective doublet. 
[0034] For AS 0 a similar formula applies: 

5 AS ° total = ZP£S°(') + AsVinit) + AS°(3 , init) + AS°(symmetry) [Eq 3] 

[0035] For Eq 2 and Eq 3 it should be noted, that the symmetry correction is zero for non-self-complementary oligo- 
nucleotides. The AH* and AS* values for Watson-Crick nearest neighbors with matched and mismatched basepairs as 
well as those for initiation and symmetry correction were taken from published data and are summarized in tables 1 
10 and 2. A standard PCR buffer contains 50 mM KCI and 1 .5 mM Mg 2 * for rapid cycling applications higher Mg concen- 
trations are often used. Mg 2 * cations have a 140 times higher stabilizing effect on duplexes compared to Na+ [13]. 
The concentrations of the monovalent Ions do not add much to this effect, at least under the conditions which apply 
for PCR [13]. Therefore the sodium equivalents are calculated as follows: 
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Na + equivalent = [Mg 2+ ] x 140 [Eq 4a] 

[0036] Alternatively a formula based on a 1 00 times higher DN A stabilizing effect of Mg 2 * cations may be used. This 
also takes into account the monovalent ions [14]: 

Na + equivalent = [Mg 2+ ] x 1 00 + monovalent ions [Eq 4b] 

[0037] The entropy is salt dependent and therefore AS 0 must be corrected if the sodium concentration is different 
from 1 M NaCI, the condition under which the parameters tabulated in table 1 were derived p]. 

AS 0 [Na] = AS°[1M Na] + 0.368 X N X In [Na*] [Eq 5] 

where N is equal to the oligonucleotide length minus 1 . 

If DMSO is included in a PCR reaction, it induces a Tm reduction of 0.6°C for each percent DMSO [15]. 

1 . Example: Calculation of Tm for a 18mer without mismatches 

[0038] As an example the Tm of the 1 8 mer 

5*-ACG ATG GCA GTA GCA TGC^3' according to Figure 4 is calculated. 
3'-TGC TAC CGT CAT CGT ACG-5' 

[0039] This example oligonucleotide has 17 n-n doublets. These are AC/TG, CG/GC, GA/CT, AT7TA, TG/AC, GG/ 
CC, GC/CG and so on, see Figure 4. 

It should be noted that the first doublet (AC/TG) is identical with GT/CA. Only the thermodynamic data for latter doublet 
the is listed in table 1 to save space. Palindromic sequences are not listed twice. From the n-n given in table 1 , the 
number of different doublets found in the above 18mer oligonucleotide are as follows: 

N AA/TT 0> N AT/TA ^» N TA/AT ^ » N CAA5T ^» N GT/CA 2 f N CT/GA N QA/CT ^ » N CG/GC ^ » N GC/CG 3> N QGVCC 1 

[0040] With data from table 1 AH° can be calculated according to Eq 2: 

E/ifAH° (0 results in 0X(-7.9)+2x(-7.2)+1 x (-7.2)+4x(-8.5)+2X(-8.4)+2x (-7.8)+1 x (-8.2H1 X 
(-10.6H3X (-9.8K1 X (-8.0)=-144.2 kcal/mol 

[0041 ] The initiation term for each direction is 0.1 kcal/mol for initiation with terminal G-C and 2.3 kcal/mol for initiation 
with terminal A-T. 

No symmetry correction is necessary since the oligonucleotide is non-self-complementary. AH° totel = -144.2 + 0.1 + 

2.3 kcal/mol = -141 .8 kcal/mol. 

[0042] AS* calculation according to Eq 3: 



5 



10 



15 



EP1 103910 A1 

Z,n,AS°(,) results in 0x(-22.2) + 2x(-20.4H1 x(-21.3>t4x(-22.7) + 2x(^2.4>f2x(.21.0) + lx 
(-222K1 X(-27.2>^3x(-24.4>f1 X(-19.9) = -3B2.2 cai/Kxmol 

No symmetry correction is necessa,y since the oligonucleotide is non-seff^omplementary. 

AS tobi = -380.9 cal/Kxmol 
[0043] Sodium equivalents calculation according to Eq 4a- 

Assuming the PGR contains 2.5 mM Mg- according to Eq 4a the sodium equh,a.ent is 2.5 m M x 140 == 350 mM = 
Erai 8 " *" according to Eq 5: 

so AS* [Na)= -380.9 cal/Kxmol + 0.368 x 1 7 x In [0.35] = -387.5 cal/Kxmol 

[0045] Calculating Tm according to Eq 1 

B^use o, the different un«s AH must be multiplied oy 1000 to caWno,. To convert from KeMn to Ceisius subtract 

25 ^™*::* ep ^ teo '^^ 



so T m = -141 800 fa/ /wo/ 

-mScal/K* moU 1.987 call K« mo lx l^M^mo/J = 627 ° C 

35 2. Example: Calculation of Tm for a 18mer with mismatches 

[0046] As an example the same 18 mer5 , -ACGATGGCARTA<-r-ATr>^o. 

mutation at the third base. Hybridization of me^d tvi^^ ™?Z? 8b0We " teken a <*>A 

wlthaGTmlsmatchinthe^(SS^Tr 8 ^ d ^ P ^ ^ * e ^ to the following duplex 
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45 



5'-ACQ ATG GCA GTA GCA TGC-3* 
3'-TGITAC CGT CAT CGT ACG-5' 



N A^0 1 N AT ^ A 2, N - AT Y > N C I CT 4 N™2 N^*^ /CG ^ 8 ° on ^ enu "*«rof different doublets are: 
so doublets affected by th^mutatlon^ N 1™ N CTA3A 1 ' * W 0. N ^ 3, N ^ 1 . The two 

pairs are equal to example 1 . CQ * r ^ ^ ° f these doub,ets o^*"* °"ce. The remaining n-n 

£21 w ^ ^/l 80 ^ teb ' e 2 AH ° can te calcu,ated ^^"3 to Eq 2: 

* P^areSdin^ 

is the same as GA/TT) AH 0 = -13 kcal/mol aou bletAH - -4.1 kcal/mol, forGAfiT (look up TT/GA in table 2, this 

!^Tom„:^ 
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No symmetry correction is necessary as the oligonucleotide is non-self-complementary. 



AH ° total = - 125 - 4 " 4 -1 - 1 -3 +0-1 +2.3 = -128.4 kcai/mol 



[0049] AS 0 calculation according to Eq 3: 

For the Watson-Crick paired bases £, n,AS0ty results in 0x(-22.2)+2x(-20.4)+1x(-21 .3)+4x(-22.7)+2x(-22.4)+2x(-21 .0) 
+0x(-22.2>+Ox(-27.2>^3x(-24.4H1x{-19.9) = -332.8 cal/K x mol. Two doublets are affected by mismatches and their 
thermodynamic properties are again found in table 2. For the CG/GT doublet AS 0 = -11.7 keal/Kxmol, for GA/TT (look 
up TT/GA in table 2, this is the same as GA/TT) AS 0 = -5.3 kcal/Kxmol. 

The initiation term is added for each direction, that is -2.8 cal/Kxmol for Initiation with terminal G-C and 4.1 cal/Kxmol 
for Initiation with terminal A-T. 

No symmetry correction is necessary as the oligonucleotide is non-self-complementary. AS° totaJ =-348.5 cai/Kxmol 
[0050] Sodium equivalents calculation according to Eq 4a: 

[0051] Conditions are the same as above: 2.5 mM Mg 2 * according to Eq 4a the sodium equivalent is 2.5 mM x 140 
= 350 mM = 0.35 M. 

[0052] Calculating AS° for a given sodium equivalent according to Eq 5: 
For a 1 8 bp duplex is N=1 7. 



[0053] Calculating Tm according to Eq 1 

Again AH* must be multiplied with 1000 to change to cal/mol. 273.15 is subtracted to convert from Kelvin to Celsius. 
The concentration of the probe is 0.1 |iM = 0.1 x 1 Qe-6 M and this must be divided by 4 because the sequence is non- 
self-complementary. 



3. Example: Calculation of ATm for a hybridization probe mutation detection assay 

[0054] In the first example a Tm of 62.7°C was calculated for a 1 8mer probe. In the second calculation a lower Tm 
of 56.2*C was calculated due to a single mismatch caused by a G->A mutation. If the former probe were to hybridize 
on the latter mutated strand, a calculated reduction in Tm of 6.5°C would result 

[0055] The thermodynamic calculations assume that the annealing occurs at pH 7.0. This is not the case in standard 
PCR buffer systems, but the resulting error seems to be small. More important is to note that mismatches in the terminal 
and penultimate position are not property reflected by the parameters in table 2 [3;10]. The thermodynamic data set 
for double mismatches (i.e. , Immediately adjacent mismatches) is still incomplete. Only for the few double mismatches 
listed in table 2 is a reliable prediction possible. The influence of the salt concentration on oligonucleotide stability is 
also still under investigation [4;5;13], Nevertheless, the methods presented here will allow a reliable estimate of hy- 
bridization probe Tm on the LightCycler but may require correction of the sodium equivalents according to the local 
PCR conditions. 

4. Validation of the n-n model on the LightCycler 

[0056] The formulas described above were incorporated into a Microsoft® Excel® spreadsheet software which is 
capable of calculating the thermodynamic data for matched and mismatched oligonucleotides under different concen- 
trations and ionic strengths. This program is available for download via http7/serveri.medikc.med.uni-goettingen.de/ 
meitcalc.htm or it may be requested by email. A total of 87 melting points resulting from eleven different published 
[16-26] and 14 unpublished hybridization probe assays where calculated and showed an excellent fit with the predic- 
tions from the n-n model, as may be seen from Figure 5. For the data shown in Figure 5 Eq. 4b was used to calculate 
the Na+ equivalents which gave a better fit when using the published Tm data from different laboratories. However, 
under the standard PCR conditions used with 2.5 mM Mg 2 * and 5% DMSO the best prediction were obtained using 



AS* [Na)= -348.5 cal/Kxmol + 0.368 x 1 7 x In [0.351 = -355.1 cal/Kxmoi 



-128400 kal I mol 



- 3 55.1 cal/Kx mol + 1.987 cal/Kx mol x 
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Eq4a. 

POST] ^notyping assays for the apolipoprotein B-3500 and factor V Leiden mutation have been performed. For 
?f mu ^ tton ? ave ^ described *at occur under the detection probe. These are the G10699A and 
££2^ 22? °' ^ ^ in B 1271 and thB G1691A - A1692C and G1689A *« *» factor V P0;»] 
EZSJ? ? f °r^ e PreVa ' enCe ° f Sin9te " UC,e0,lde P°*™PN™ the human genome [29;30] raises the 

?°2 ^ P 088 "" 8 mUtati0nS Under a 19mer probe led us that no mutation under a 

probe binding to the same extent [12J. These mutations will not be discriminated from each other. Fig 3 shows how 

JZJSZT J?™" 1 ^ T tem ^ deStabili2ed by S Si "9 te base mute tion. The most |qj£ SI Z 
probes with higher GC content show less Tm shift due to a single mutation. Another principal point is that shorte 
probes showhigherTmshitts because the destabilizing impact of the am^antoJLJL^JmwTZ 
duplexes. The them,odynamic n-n model provides the theoretical basis for the stabiltty of matched and Sri 
snfw^Sa^r 80 " ° f eXpenmenta ^ derived P— «■» thosTcafculated ST^nTrS 
Using adetection probe that has been proper* optimized according to mis mc^ 

of every point mutation should be possible. However, it is possible that different mutation?Ler a probe res^S" 
same Tm reduction. Important factors for the design of hybridization probes are: (I.) Absence of sel - or 
cross dimertzatjon of probes and amplification primers. (II.) Equal meiting points o! detection and a^^fand 
absence of false pnming sites on the amplicon for hybridization probes used for quantification. (Ill ) Use of a detection 
probe wrth a maximum Tm shift caused by the Investigated mutation (see below) and an anS w^T^m 
approximately 10% higher than that of the detection probe for mutation detection. Provided ft is noTin fte^eTlt^Le 
or ultimate position the place where the mutation is located under the probe is not important The "JmEZ meS 
pointshiftwift be affected by me type of mutation and the neighboring bases but not SS^^S^lSSS 
The commonest mutations are G->A mutations [31). However, the resulting G:T mismatch underlne fSTJSL 

Se rl^L 9 T" 3 * ° 0BHr t098ther Wfth G:A and G:G mfematohes » « * therefore adSZleto^Z 

Sf^St P °1 E !. antiSen8e Strand ThB reSUlting ^P^tary C:A mismatch has a higher destabC 

SSSSHHn * ?! t8Cted i ,V - ) Th8 93P betWee " the fluoro P^es should be as small as possible since the 

k inVerSely 10 th9 SiXth P° Wer 0f me distance betwee " the donor and the acceptor dye 

[32J. One to five basepairs gap has been reported to give good results [1]. 
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Table 1. 



Unified parameters for oligonucleotide AH* and AS" in 1 M NaCI [5]. " ~ 1 

(Nearest neighbor dimer duplexes are presented with a slash separating strands in antiparaliel orientation e q AT/ 
TA means 5'-AT-3' Watson-Crick paired with 3'-TA-5\ ) 


Sequence 


AH* 


AS* 




kcal/mol 


cal/Kxmol 


AA/I | 


-7.9 


-22.2 


AT/FA 


-7.2 


-20.4 


TA/AT 


-7.2 


-21.3 


CA/GT 


-8.5 


-22.7 


GT/CA 


-8.4 


-22.4 


CT/GA ~ ' "~ 


-7.8 


-21.0 


GA/CT 


-8.2 


-22.2 


CG/GC 


-10.6 


-27.2 


GC/CG 


-9.8 


-24.4 


GG/CC 


-8.0 


-19.9 


Initiation with terminal G-C 


0.1 


-2.8 


Initiation with terminal A-T ~~ 


2.3 


4.1 


Symmetry correction for self complimentary oligos 




-1.4 



Table 2. 



Parameters for mismatched oligonucleotide AH U and AS* in 1 M NaCI " " n 

\n»Z^ ^ 9 ^ ^ NearGSt neighb ° r dlmer dUp,exes are ^ a s,as " separating strands 
SK? .'S^ e ' 9 * AT/FG mismatehed wHh 3--TG-5V Mismatches are in^e second pc! 
srtion of the listed doublet except for double mismatches. These parameters are not valid for the calculation of 
mismatches in the terminal or penuftimate position.) 


Sequence 


AH° 


AS* 




kcal/mol 


cal/Kxmol 


Internal G-A mismatches [9J 






AA/TG 


-0.6 


-2.3 


AG/TA 


-0.7 


L "2.3 


OfiJGG 


-0.7 


-2.3 


CG/GA 


-4.0 


-13.2 


GA/CG 


-0.6 


-1.0 


GG/CA 


0.5 


3.2 


TA/AG 


0.7 


0.7 


TG/AA 


3.0 


7.4 


Internal C-T mismatches [8] 






AC/TT 


0.7 


0.2 


_ AT/TC 


-1.2 


-6.2 


CC/GT 


-0.8 


-4.5 


CT/GC 


-1.5 


-6.1 
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Table 2. (continued) 





Sequence 


AH° 

kcal/mol 


AS* 

cal/Kxmol 


5 


Internal C-T mismatches [8} 








GC/CT 


2.3 


5.4 




GT/CC 


52 


j 13.5 


10 


TC/AT 


1.2 


0.7 




TT/AC 


1.0 


0.7 




Internal A-C mismatches [7] 








AA/TC 


2.3 


4.6 


15 


AC/TA 


5.3 


14.6 




CA/GC 


j 1.9 


3.7 




CC/GA 


0.6 


-0.6 


20 


GA/CC 


5.2 


14.2 




GC/CA 


-0.7 


-3.8 




TA/AC 


3.4 


8.0 




TC/AA 


7.6 


205 


25 


Internal G-T mismatches [3] 








AG/TT 


1.0 


0.9 




AT7TG 


-2.5 


-8.3 


30 


CG/GT 


-4.1 


-11 ;7 I 




CT/GG 


-2.8 


-8.0 




GG/CT 


3.3 


10.4 




GT/CG 


-4.4 


-12.3 


TG/AT 


-0.1 


-1.7 




TT/AG 


-1.3 


-5.3 




Internal A-A mismatches [10J 






40 


AA/TA 


\2 


1.7 




CA/GA | 


-0.9 


-4.2 




GA/CA I 


-2.9 


-9.8 


45 


TA/AA 


4.7 | 


12.9 


Internal C-C mismatches [10] 








AC/TC 


0.0 


-4.4 ! 




CC/GC 


-1.5 


-7.2 


50 


GC/CC 


3.6 


8.9 




TC/AC 


D.l 


16.4 




Internal G-G mismatches [10] 






55 


AG/TG 


-3.1 


-9.5 


CG/GG 


-4.9 


-15.3 




GG/CG 


-6.0 


-15.8 
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Table 2. (continued) 


Sequence 

Internal fi-f^ mkmatrhao msyi 
iiiu;iiiai \j v3 1 1 iio rriaicn es |1U| 


Lr > al/mrvl 
KUUI/mOl 


AS* 

cal/Kxmol 








1.6 


J 3.6 








AT/TT 


-2.7 


-10.8 


CT/GT 


-6.0 


-15.8 


« I/O 1 


-2.2 


-8.4 


1 l/AI 

Double mismatches [3] 


05 


-1.5 


GG/TT 


5.8 


16.3 


GT/TG 


4.1 


9.5 


TG/GT 


-1.4 


-6.2 



Claims 

1 - A method of automatically selecting oligonucleotide hybridization probes for 
match in a DNA duplex, the method comprising the steps of: 

- defining the mutation to be detected by 



detecting a mutation causing a mis- 



a DNA sequence of a mutant DNA strand and 

a DNA sequence of a wild-type DNA strand corresponding to the mutant DNA strand; 
m^ 

defining a range from ml to nl of possible values of the variable Lt • 

TrH!Z an I" 1696 ' Va ?^! e 12 f ° r 8 SeC ° nd ' ength ° f *° onsonuoleotide hybridization probes from the place 
of mutation in a second direction opposite to the first direction- P 
defining a range from m2 to n2 of possible values of the variable L2 
selecting one or more kind(s) of the oligonucleotide hybridization probes from 

probes which have an one-to-one hybridization fit to the mutant DNA strand 

DNA^, h haVea " ° ne ~ t0 "° ne hybrWiZati0n mt ° ^^"^DNAstmnd'correspondingtothemutant 

probes which have an one-to-one hybridization fit to an opposite sense mutant DNA strand which fe th» 
complementary DNA strand to the mutant DNA strand, and the 
probes which have an one-to-one hybridization fit to an opposite sense wild-type DNA strand whw, I* th« 
complementary DNA strand to the wild-type DNA strand corresponding to themZS^ sC 

?s?oS^ 

def, K"l 9 f t he ™ od y namic nearest-neighbor model for calculating a melting point of a certain hvbridtertion 
probe hybndizedtoaDNAstrand.themelting point being a temperature at vvhi^ 

calculating, on bas.s of the thermodynamte nearest-neighbor mode., for ail possible ^M^t*^ 

a first melting point of the respective probe hybridized with its complementary mutant DNA strand or oo 
posite sense mutant DNA strand, respectively, ° 0r op " 

a second melting point of the respective probe hybridized with its complementary wild-type DNA strand 
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or opposite sense wild-type DNA strand, respectively; and 

a temperature difference between the second melting point and the first melting point; and 

- ranking the possible hybridization probes with regard to the temperature difference. 
The method of claim 1 , and further comprising the steps of: 

- defining a temperature range of a hybridization assay to be employed for detecting the mutation using the 
oligonucleotide hybridization probes; and 

- checking for each of the possible hybridization probes whether both the first and the second melting points 
are within the temperature range of the hybridization assay; wherein the step of ranking the possible hybridi- 
zation probes with regard to the temperature difference includes: 

- ranking only those possible hybridization probes which have both their melting points within the temperature 
range of the hybridization assay. 

The method of claim 2, wherein the step of defining the temperature range of the hybridization assay to be employed 
for detecting the mutation Includes: 

defining the temperature range within a range from 20 to 80 °C. 

The method according to any of the claims 1 to 3, and further comprising the step of: 

- defining at least one physico-chemical parameter of a hybridization assay to be employed for detecting the 
mutation using the oligonucleotide hybridization probes; wherein the step of defining a thermodynamic nearest- 
neighbor model includes: 

defining the thermodynamic nearest-neighbor model dependency of the at least one physico-chemical param- 
eter of the hybridization assay. 

The method of claim 4, wherein the step of defining at least one physico-chemical parameter of the hybridization 
assay to be employed for detecting the mutation includes: 

defining at least one ion concentration as one of the physico-chemical parameters. 

The method of claim 5, wherein the step of defining at least one ion concentration as one of the physico-chemical 
parameters includes: 

- defining a Na+ cation equivalent concentration. 

The method according to any of the claims 1 to 6, and further comprising the steps of: 

- defining an integer number N of the hybridization probes to be automatically selected; and 
making a list of the N top ranked possible hybridization probes. 

The method according to one of the claims 1 to 7, wherein the step of ranking the possible hybridization probes 
with regard to the temperature difference includes: 

- ranking the possible hybridization probes separately for each of the kinds of the oligonucleotide hybridization 
probes. 

The method according to one of the claims 1 to 8, and further comprising the step of: 

- checking, for each of the possible hybridization probes, whether it is able to hybridize either to the mutant DNA 
strand and the opposite sense mutant DNA strand, or to the wild-type DNA strand and the opposite sense 
wild-type DNA strand in more than one way; 

wherein the step of ranking the possible hybridization probes with regard to the temperature difference includes: 

- ranking only those possible hybridization probes which are neither able to hybridize to the mutant DNA strand 
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10. The method according to one of the claims 1 to 9, and further comprising the steps of: 

" 2£S£ hybrid ' Zati0n ^ ^ " fe ^ t0 thB ™**>" * *e 

wherein the step of ranking the possibie hybridization probes with regard to the temperature difference Includes: 
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Co cn 

> -N 



GT/CA 



CA/GT 



CT/GA 



CT/GA 



CA/GT 

Inititiation with 
terminal A-T 



1 


| 
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O 
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►J* 






i T 






u J 
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Q 
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Fig. 4 
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